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In this second paper of tlie series we calculate the stress tensor of excited matter, created by 
"debris" of high energy collisions at the boundary. We found that massive objects ( "stones" ) falling 
into the AdS center produce gravitational disturbance which however has zero stress tensor at the 
boundary The falling open strings, connected to receeding charges, produce a nonzero stress tensor 
which we found analytically from time-dependent linearized Einstein equations in the bulk. It 
corresponds to exploding non-equilibrium matter: we discuss its behavior in some detail, including 
its internal energy density in a comoving frame and the "freezeout surfaces" . We then discuss what 
happens for the ensemble of strings. 



PACS numbers: 



I. INTRODUCTION 



This is the second paper of the series, and the first one 
[l| (to be referred to as [I] below) had rather extensive in- 
troduction. Therefore we just briefly reiterate here what 
are the main goals of this study. 

Holographic description of A/" —4 SYM theory in strong 
coupling regime can be achieved via AdS/CFT corre- 
spondence, which relates it to string theory in AdS^ x 
space, in classical supergravity regime. Large number of 
applications use this tool to study properties of strongly 
coupled Quark-Gluon plasma: but most of them are done 
in a static setting, with fixed temperature via Witten's 
AdS black hole metric. 

High energy hadronic collisions in QCD are very dif- 
ficult problems. They are time-dependent and include 
non-equilibrium physics (only collisions of heavy ions can 
be approximated by hydrodynamics and locally equili- 
brated QGP). On top of that they involve different scales 
and different coupling regimes. Pure phenomenological 
approaches were developed long ago, such as e.g. Lund 
model Q which are based on a picture of QCD strings 
stretched by departing partons. More recent approach 
- known as the Color Glass picture - was proposed by 
McLerran and Venugopalan [2] who argued that since 
fluctuations high energy collisions lead to large local color 
charges (in the transverse plane) they would lead to pro- 
duction of strong color fields. Those are treated by clas- 
sical Yang-Mills eqn in a weak coupling regime. 

Arguments suggested recently put forward a view that 
QCD has certain "strong coupling window" . In partic- 
ular, Brodsky and Teramond j4j] have argued that the 
power scaling observed for large number of exclusive pro- 
cesses is not due to perturbative QCD (as suggested orig- 
inally in 1970's) but to a strong coupling regime in which 
the running is absent and quasi-comformal regime sets in. 
Polchinski and Strassler [sj have shown that in spite of ex- 
ponential string amplitudes one does get power laws scal- 
ing for exclusive processes, due to convolution (integra- 
tion over the z variable) with the power tails of hadronic 



wave functions. One of us proposed a scenario [6[ for 
AdS/QCD in which there are two domains, with weak 
and strong coupling. The gauge coupling rapidly rises 
at the "domain wall" associated with instantons. Pion 
diffractive dissociation is a process where a switch and 
weak coupling domains are observed: and cross section 
behavior is consistent with Polchinski-Strassler expres- 
sion and expected coupling change. Last but not least, 
such approach looks now natural in comparison to what 
happens in heavy ion/finite T QCD, where we do know 
now that at comparable parton densities the system in- 
deed is in a strong coupling regime. 

Accepting the Color Glass picture as an asymptotic 
for very high parton density and large saturation scale 
Qs — > oo, one should ask what should happen in the 
case of saturation scale falling to intermediate momenta 
Qs ~ .3— 1.5 Gey associated with "strongly coupled win- 
dow" . This is the issue we address in this work, using the 
AdS / CFT correspondence in its time-dependent version, 
as a tool to describe the evolution of the system. 

The setting has been discussed in [I] , where we exten- 
sively studied how exactly the "debris" produced in a 
collision - particles or open strings - are falling under 
gravity force into IR (the AdS throat). In this work we 
do the next technical step and calculate the back reaction 
of gravity by solving the linearized Einstein eqns for met- 
ric pertubations, deducing the space-time dependence of 
the stress tensor Tp,y(x) of excited matter observed on 
the boundary. 

The general setting is in fact rather similar to the Lund 
model: except that strings are departing from our world 
{z = boundary) rather than breaking. Technically our 
work is a development along a line actively pursued by 
many authors. (Early work in a different scenario focused 
on the effective stress tensor on a branej?!)- In partic- 
ular, it can be considered a continuation of our recent 
work [s^l in which we calculated static (time-independent) 
stress tensor associated with Maldacena's static dipole. 
It is also similar to recent AdS / CFT calculations of a hy- 
drodynamical "conical flow" from quenching jet in QGP 
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The process we describe resembles what happens in 
heavy ion coUisions, but with very important distiction. 
In the setting of this paper we treat "debris" as smaU per- 
tubation, solving linearized Einstein eqns in pure AdS^ 
background. Therefore there is no black hole and/or tem- 
perature in this work, and our "mini-explosion" produce 
matter which is not equilibrated and the resulting stress 
tensor cannot be parameterized hydrodynamically. 

To get all that one needs to proceed to nonlinear grav- 
ity and a non-linear process of black hole formation: the 
problem which we hope to attack elsewhere. 



II. SOLVING LINEARIZED EINSTEIN 
EQUATION 

We want to solve linearized Einstein equation in AdS^ 
background in Poincare coordinates, with standard back- 
ground metrics 



, o -dt'^ + dx^ + dz'^ 
ds-^ = ^ 



tained via variation of their action over the metric 
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here g^v and gind.afS denote the AdS metric and the in- 
duced metric on the string worldsheet, respectively. 

In the next section we find an expression for Green's 
function, which will provide hmn for any given source 
Smn- We will then exctract an expression for the coeffi- 
cient of the term in Taylor series of hmn at the bound- 
ary, which by the rules of AdS/CFT correspondence gives 
us the boundary stress tensor. 



An axial gauge is chosen for metric perturbation, so 

/izM =0 

The linearized Einstein equation are well known and 
the tactics used in its solution are discussed in Q: the 
present case is only different by appearence of time 
derivatives. We put it into the form 



1 z 
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where D = z'^ {—^f + d'i + d^), the indices are 0-3 and 
the r.h.s. is the generalized source 



JQ 

containing not only the stress tensor of the source 



T 



(3) 
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but also the following combinations of perturbation met- 
ric which can be easily found from the eqns: 
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The source term for different objects is standard, ob- 



III. THE GREEN'S FUNCTION FOR THE 
LINEARIZED GRAVITY IN AdSs 

The Green's function we need satisfies the following 
eqn: 

2 

^{dl ~ + a|)G(z, z') - 2G(z, z') + ^dzG{z, z') 



6{z - z')S{t - t')d^-^\x - x') 



(7) 



and the solution to ([2]) is then given by hmn{z) = 
J G{z, z')smn{z')dz' . Thus G{z,z') should satisfy the 
same boundary condition as h{z). Fourier transforming 
4-dim part of ([7]), we have z- dependent eqn 



y + _ fc2)^(^^ _ 2G(z, z') + |5,G(z, z') 



5{z-z') 
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where G'=(z,z') = / G{z, z')e-^'^*+'''^dtd^x 

([U can be solved in terms of Bessel functions: For 
\uj\ > /c, the solution is a linear combination of Bessel 
functions of the first and second kind. We impose the fol- 
lowing boundary condition: at z = 0, G{z, z') = 0(h{z) = 
0), at 2; — > oo, G{z, z'){h{z)) contains outgoing wave only, 
i.e. the wave propagates from the source to infinity ^0|. 
The solution is composed of two homogeneous solutions: 



'^(Z,Z ) - i 7 



Z < Z' 



B{J2{\z) - isgn{uj)Y2{Xz)) z > z' 



(9) 



with A = \/lo'^ — k'^, A, B is fixed by matching the func- 
tion itself and its first derivative at z = z': 
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J2(Az') 



(10) 



For k > \uj\, the solution can be built from Modified 
Bessel function. We choose the no exponential growth 
boundary condition at z ^ ooQ. The solution is given 
by: 



G(z,z') = { 



ChiXz) z <z' 
DK2{\z) z> z' 



(11) 



with A — \/k? 



^ _ -2_ff2(Az') 

^ ^ -2l2{\z') 



(12) 



It turns out the solution can be organized in a compact 
form using properties of Bessel hmction: 

-|r/2(«Az<)iir2(«Az>) w > 0, |w| > fc 
G{z, z') = { -jrl2{-i\z<)K2{-i\z>) w < 0, |w| > A(13) 

-^h{Xz^)K2{\Zy) \U\ < k 

where z< = min{z, z') z> = max{z^ z') 

Doing the inverse Fourier transform: 
(2^)<i / {z, z')e^'^^'^^^ du}(Pk, we obtain a retarded 
propagator for the metric: (See Appendix A for the 
evaluation of the integral, a retarded scalar propagator 
was found in [T^ . Other scalar and graviton propagators 
were also found in 13] with slightly different boundary 
conditions) 



12iz 
(2^ 
1 



1 



(^2 _ ^2 



(t2 _ J.2 _ ^2 + 

e{t-r) 



(14) 



Several comments of the propagator are in order: 
(i)The theta function implies the propagator acts inside 



the lightconc t 



> and moreover is re- 



tarded t > r > 0, which we indicate by the subscript 
R. It is also consistent with the outgoing boundary con- 
dition. Note that the propagator is also Lorentz invari- 
ant. (ii)The propagator relates the coefficient of met- 
ric perturbation Qmn and the source Smn in the following 
way (assuming Smn{z) does not contain z° and terms): 

Qmn{t',x') = J PR{t' - t,x' - X,z)Smn{z,t,x)dzdtd^X 

The primed and unprimed coordinates correspond to 
boundary and bulk respectively. (iii)The propagator is 
dynamical. For static source, one can perform the t- 
integral to obtain a static propagator, which agrees with 
the one obtained in [8]. 



IV. THE STRESS TENSOR OF A FALLING 
STONE IS ZERO! 

First we want to study the stress tensor induced by a 
falling stone. The tra jectory of stone in AdS background 
is: z{t) = y^z^+t"^ = z ioT t > 0. Plugging the trajec- 
tory, it is not difficult to obtain the generalized source: 



Hi'^mS'^^^x) 



3z^ 
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(z2 + 2t^)e{z - z) 



(15) 



Convolute the generalized source with the propagator 
Pr- We use the substitution: dx^ = — d = 9 x'^ 
(where Xm = tiX^) to simplify the calculation. We use 
partial integration in the first identity so that the deriva- 
tive acts on the propagator, which we indicate by an 
over left arrow. The second identity is due to the fact 
Pr = Pr{x'„^ — Xm). After that we can see that the a;* 
integral kills the S^^^{x) and the z-integral involving ei- 
ther delta function or theta function can also be done 
easily. Finally we are left with the i-integral: 
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2{t' t) ^x'\ 2,p it'- if -it' - ty:, ^ 

)^^tP+[_^t'-t)x', x'^r IX 



(z^ + 2t^)P 
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with 



P 
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(2^ 



1 



((i'-f)2-r'2-z2+ie)4 
1 



((<' - <)2 - r'2 - z2 - iey 



e{t' - 1) 



(17) 



The i-integral is evaluated via residue theorem. No- 

tice P contains a fourth order pole at t = 2p — ~ i^i 

the denominator, only the coeffcients of and i"* in the 
overall numerator are relevant for the integral. Summing 
up the matrices, we end up with a vanishing result. [2l|. 
It is probably not a coincidence that trF'^ 12] calculated 
from dilaton perturbation vanished as well. Both calcula- 
tions are in linearized approximation and correspond to a 
falling stone. Nevertheless the nonlinear version of stress 
tensor induced by the same object obtained in [3, HB] 
seems to contradict our result, which deserves further 
study. 



4 



V. THE STRESS TENSOR OF A FALLING 
OPEN STRING 

We want to study the stress tensor by a falling string. 
A scaling string profile is obtained in [I] for a separating 
quark- antiquark pair, provided the separating velocity 
is not too large: v < 0.6. We briefly recall the scaling so- 
lution. The quark(antiquark) moves along the trajectory 
X — ±vt. The string profile is given by: 




T and y are proper time and space-time rapidity. The 
y = Y limit of lfT8|) relates the parameter /o and the 
quark rapidity Y = arctanh(v). It is also very useful to 
write down the EOM of f{y): 



/' = 



^V{V-E^) 
E 



(19) 



withy = /4-/2,/4-/2_i^2^0. 

We want the source term due to the scaling string. It 
is convenient to switch to a different parametrization: 



z — z(t,x) , = 



(20) 



where x and x± represent longitudinal and transverse co- 
ordinates respectively. The above parametrization leads 
directly to 



1 k'z 
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2(1 - z 



4)) 



(21) 



(In the matrices here and below we only show the 
nonzero enties: the adopted order of coordinate indices 
is ^^a:^a:^x3. ) With the help of string EOM, (15) 
of llj, h can be expressed in a very compact form: 

h = -1 2^<5 (x^) S (x^) ^(z^ - z^)0{z - z) We 

also record the generalized source for later reference: 
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(22) 
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With the source now at hand, we proceed to the cal- exphcitly, we find the z-integral and xj^-integral can be 
culationjsf stress tensor. ^We use similar substitutions as done easily. We arrive at the following result: 
before: dx = ~dx — d'^. Performing the derivative 



with 
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which is just the integrated propagator. The four ma- with 
trices in the expression above we will refer to later as 
I, II, III, IV, respectively. 

Here we replaced the theta function of Pr by 6{t' — t). 
It is justified since the ±ie prescription encodes deriva- 
tives of the delta function, and the theta function picks T'ch{y' — y) ± /T''^ch'^{y' — y) — (t'2 — a;'2)(l \) 

up only one pole of the propagator. — jL '. ~f26) 

In order to plugin the scaling solution for the string, ^ — yi 

we return to t, y coordinates: 



Zt 



chy shyf 

f P 
shy _ chyf 

f P 
rdrdy 



dtdx 



The source together with the integration measure has 
one of the following simple r-dependence: t,P,t^. The 
propagator now is 



P 
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This propagator as a function of r contains two fourth 
order poles, so the r-integral is calculated by the residue 
theorem. Note that the theta function picks up only one 
pole at r = T_. Since our r-integral extends from zero 
to infinity, we must have a positive t_ in order to have a 
nonvanishing result, which requires r'2 — a;^ = t'^ — r'2 > 
0, precisely the condition that the observer must stay 
inside the lightcone. Since the quark and the antiquark 
are emerging from the space-time origin, the stress tensor 
is indeed expected to vanish outside the lightcone. 



(27r)2 ((1 
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Completed the r integral and replacing ^ by , 
we convert Qmn to the boundary stress tensor Tmn (com- 
(25) pare @). We thus get our final result 
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A. A field near one charge 

We first consider the stress tensor close to the quark. 
Note the quark moves with velocity v, the observer 
should also move in order to stay close to the quark. 
Thus it is convenient to switch to rest coordinate of the 
moving quark, yet we still stay in the rest frame of the 
system. The rest coordinates of the quark, indicated by 
a tilde, relates the original coordinates in the following 
way: 

i = chYt' - shYx' 
i ^ chYx' - shYt' 



X2 

X3 



Xn 



(29) 



We expect to obtain the field by the quark only, pro- 
vided we are sufficient close to the quark, which corre- 
sponds to the limit f 3> r , f = ^Jx^ + x\+ x^ 0. 
Comparing the string profile for the stretching dipole 
with that for a single quark, we claim the leading order 
field near the quark receives contribution from the quark 



end of the string, which corresponds to the integration of 
y near y = Y in (p7)) . If we instead do the integral in / 
with dy — jr, we may only focus on large / integration. 

Note f ^ ^ y — Y ^ -p. To the leading order, we 
can simplfy replace y by F, which leads to the following 
relations: 
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Performing the / integral near 00, and take the limit 
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i ^ f ,r 0, which is essentially a small f expansion, 
we find the leading order field given by source I and IV 
diverges as 0{^), while source II and III only yield sub- 
leading contribution 0( Jj). We display the LO field near 
the quark as follows: 



the first one is source II and III convoluted with P4, the 
second correspond to source IV and P4, and the third one 
comes from the convolusion of source I and P5 . Collecting 
all of them, we find the following result: 
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JXX3 
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12f6 



(31) 



with 7 = chY , 7/3 = s/iK . 

This does not look very nice at first glance, actually 
it is just the stress tensor of a static quark boosted to 
a frame moving with velocity —v. It is clearly traceless 
and divergence-free. We confirm that the LO field near 
the quark contains contribution from the quark only. 

Next we would like to extend the result to NLO to 
include the effect of the antiquark. Note there are two 
possible corrections relevant for NLO: correction to the 



72/3252 

-72/3x2 

-7/3x^2 
-7/35x3 

/ 272/325 
272/35 
-7/3^2 
-7/3^3 



-72/352 

2-2 

7'^x 

7552 

IXX3 

-272/35 
2725 

7^2 
7^3 




(33) 



source Ay = y — Y 
shYAy, shy =^ shY -\ 
the propagator: 



thus chy = chY 
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Al _ E 
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chYAy. The other is correction to 



±6 



We had found with satisfaction that this result is indeed 
traceless and divergence- free (to order 0{jt)). 

After this result is boosted to a frame moving with 
velocity +v, it reproduces the NLO field near the quark of 
a static dipole with the identification -p- ^ p- « ^^{L 
is the quark- antiquark seperation). It confirms that close 
to one charge it is not important to this accuracy what 
the other charge is doing. In fact in the quasi-static limit 

E 



seperation at time t. 



Ip-, where L is the quark- antiquark 
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12i 
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((1 - 72 )t2 + t'2 - 2TT'ch{y' - y) - x'l ± itf ' 3/ 
■■■]6{r'-T) 

[Pi + P5 - 



B. The slow-moving limit 

Since the stretching string solution depends on only 
one parameter w, the ends velocity, which is bounded 
from above by the critical velocity. One interesting limit 
in which calculations can be pushed a step further is 
w — > 0, which correspond to slow motion. In practice, 
it corresponds to expantion of the stress tensor in inverse 
powers of /q. 

We start with considering the large /o limit of (|27p . 
Define r/ = 77 ^ 1 such that the range of ij is inde- 
pendent from /q. The large fo expansion of y is a little 
complicated: 
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±1 



Y) 
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After relatively lengthy calculation and comparison we 
find that the NLO correction is composed of three pieces: 
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where G{ri) = -^<^^'^H2g(^^'^) 

With the asymptotic expansion of y and /o, we are 
ready to proceed to the stress tensor. It seems at first 
glance the leading order is of 0(-^), given by source IV. 
Actually the prefactor, which is an integral of r\ vanishes. 
The order 0{j^) does not contribute either due to the 



symmetry y ^ —y. Finally we have to extend the calcu- 
lation to order O(-i-). Expand all the relevant quantity 

Jo 

in /o, and keep the order 0(73-) in the result of stress 

JO 

tensor. It is a quite lengthy but straight forward calcula- 
tion. The result is displayed as follows(we have omitted 
the prime in boundary coordinates): 



,41 



Ttt = [(10a + 5ei + 5e2)r^t^ + 45eia;^r^ - 35eii^a;^ - (9e2 + 2/ + 6a - 24c + 7ei)r^] 

/o ^ 

2\/X 2x 

Ttx = 75-^-9 [(45ei + 1562 - 90d - 30c)rH'^ + ISeixV^ - 105eit'^x^ - (Sea + 7ei + 2/ - 6c - 18d)r' ^ 
Jo''' 

Tte^ = .?^^[(15e2 + 15ei - 30c)rH'^ - lOSeii^a;^ + 15eix^r^] 
JqT^ f 

T^x = -7^^ [(20a - 306 + lOei - 60rf)r*i^ + (-12a + 186 - 6ei + 36rf)r^ + (420rf - 175ei - 2,f)e2)rH'^ 
-|-(-180rf + 1562 + 6561 + 10/)a;V^ + 315eixH'^ - 1056ia;V] 

Txx, = If"^ ^"fi""' [(42d - 21ei - 7e2)fr' - 21e,x\^ + 63e,tV] 
Jo''' 

Tx,x, = |$^(5«i' - Sar')Sij - [(7ei + 7e2)r2i^ - (362 + ei + 2/)r^ + 21e^x'r' 

-eSeixH'^] (35) 



with 



x^ "H X2 ~\~ x^ 



a= / dn = 0.5991 

f°° 1 
b= ^^dr? = 0.3594 

1 V 7?V»? - 1 



drj = 0.4493 



_ r G(7?) 



drj = 0.0899 



61 = / / f Ginfdv = -0.1797 
3-rf 

2v^;^43T(^2 + l)^6' 



62 = ^ / r^^ drj = 0.4793 
/ = ^ ^ "'^ ' " ^ dr] = 0.7189 



(36) 



r 



Several comments about the result are in order: (i)the lightcone. The discontinuity of the stress tensor on the 
result applies for arbitrary point on the boundary, lightcone is a consequence of the discontinuity in source 
i.e. general t, x, X2,X3, provided the point lies inside the at t = (ii)trace and divergence of the stress tensor van- 
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ish for any points away from the trajectory of the dipole 
ends, which is of course imphcitly assumed in our calcu- 
lation. (iii)If we consider the limit f ^ r, which amounts 
to keeping only the highest power of t. Recalling the 
quasi-static limit: jj ~ (5^)^ = (o)'^' where L is the 
dipole size at time t. While for the case of static dipole: 
= (r^)'^' fi'^d tl^s stretching dipole result agrees 
with static dipole in the double limits: v ^ 0,t ^ r. The 
numerical factors matches as well, (vi)the agreement of 
quasi-static result and NLO near field with those of static 
dipole by the simple identification: L = 2vt seems to 
suggest that the vacuum-quark(antiquark) interaction is 
instantaneous. 

We plot the energy density T'^'^ and the momentum 
density (energy flow) as a function of the spatial co- 
ordinates at three different times in FigH] and Fig[21 re- 
spectively. We observe that although the shape of energy 
distribution becomes more alongated with time, reflect- 
ing changing shape of the string, the shape of the mo- 
mentum flows seems to stay the same, with an interest- 
ing "eight" shape or forward-backward depletion. Small 
arrows display the direction of the energy flow. Although 
they overall show outgoing explosion away from the ori- 
gin (the collision point), one can also see some "cumula- 
tive" flow with jets converging along the collision axes. 

VI. THE ENERGY DENSITY OF MATTER IN 
COMOVING FRAME AND FREEZEOUT 

It is illuminating to ask if the stress tensor we obtained 
can or cannot be described by some hydrodynamical flow. 
The latter is widely used in describing heavy ion collisions 
[13, [III • More precisely, the question is if our stress tensor 
([55]) is that of a flowing ideal liquid 

Ti_,i, = (e + p)u^u^ + pg^^ (37) 

where is the 4- velocity of the liquid and e and p are the 
energy density and the pressure. (Tracelessness would of 
course demand that e — 3p) It is not difficult to show 
that this is not the case: the structure of our answer is 
richer than this simple form. 

Nevertheless, it is still possible to define a "comoving 
frame" of matter at any point, in which the (boosted) mo- 
mentum density Tq^ vanishes. The (boosted) local value 
Tqq component is the energy density in such a comoving 
frame, which we denote by e like for a liquid. We use the 
contour of e to define the freezeout surface. 

However the (boosted) local values of other spatial 
componets are in general unrelated and can be viewed as 
"anistropic pressure" of non-equilibrium matter. Need- 
less to say, it remains unknown what combinations of 
fundamental fields of the N—A theory - the gauge one, 
the fermion or the scalars - participate in this fiow of pro- 
duced matter: to learn that one should do "holography" 
for many more operators on the boundary. 

Recall that the stress tensor in different frames are 
related by T^^ = S^S'^T,,^ where 5^' = are Lorentz 




t=10r 





FIG. 1: (color online) The contours of energy density T"", in 

unit of 3 2 , in X1—X2 plane at different time. The three plots 

are made for t = r, t — lOr and t — 50r from top to bottom. 
Note the quark/ant iquark is at xi — ±vt. In the slow moving 
limit, they are nearly at the origin. The magnitude of T"" 
is represented by the color, with darker color corresponding 
to greater magnitude. As time increases, the shape of the 
contours gets elongated along xi axis 



boost matrix. The primed quantities correspond to new 
frame. Therefore the aim is to find such boost which kills 
all the (0, i) components. 

In practice, it is achieved numerically by the following 
recipe: 

We pick up any point inside the lightcone, calculate the 
eigenvalue and associate eigenvectors of the correspond- 
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fixed t 



10- 



p20 



.0.25 



xl 



0.5 



FIG. 3: (color online)the profile of e, in unit of -Itt', at t = 1 

Jo 

with r ~ 0.2 — 1. The evolution of the shape of contour, i.e. 
the freezeout surface is contained in this plot. The contours 
with large r (small e) are nearly spherical while the contours 
with small r(large e) are elongated along xi axis 



setting, this translates to the fohowing: At early time 
the local energy density contour is nearly spherical and 
at late time it gets elongated along xi axis. 



FIG. 2: (color online) The contours of momentum density T"', 
in unit of / \ , in xi — X2 plane at different time. The three 

Jo ^ 

plots are made for t = r, t = lOr and t — 50r from top to 
bottom. Note the quark/antiquark is at x\ — ±vt. In the slow 
moving limit, they are nearly at the origin. The magnitude 
is represented by color, with darker color corresponding to 
greater magnitude. The direction of the momentum density 
is indicated by normalized arrows 



ing stress tensor matrix. Out of the four eigenvalues, 
one is selected to be the local energy density based on 
its eigenvector (See Appendix B for a short explanation). 
Figl3] is a surface plot of e profile in spatial coordinate. 
The plot is made for t = 1. It shows a nearly spherical 
shape for contour with large r and an elongated shape for 
contour with small r. By virtue of conformality of the 



VII. THE STRESS TENSOR OF MULTIPLE 
STRINGS 



A simple extension of what is done above is to con- 
sider many coUidng quark-antiquark pairs uniformly dis- 
tributed in the transverse plane. Every pair with the 
same transverse coordinate is connected with a string. 
Let's assume quark and antiquark only interact pairwise, 
which in the dual picture means the strings do not inter- 
act with each other. As a result, the overall stress tensor 
induced by the mutiple strings, in the linearized approxi- 
mation simply amounts to integrating p5p over the trans- 
verse coordinates. Note in order to preserve causality, the 
integral is done for < x\ = r'^ — x'^ Kt^ — x^ We display 
the result as follows (we omit a factor of transverse string 
density, which does not alter space-time dependence of 
the stress tensor): 
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A 2 1 



_ 8^/A. 
20 2 

yei - -/ + 8c- 3e2 



2a) 



t 



-3 +ie2 + 2a-4ei)-3J 



8\/A r „ 2;3 -20 , _ , 2 , „ , x 
73^ [I2ei- + (— ei + 4c + 12d + -/ - 2^2)^ ^ 



2 2 1 t 

(-61 + 2c- -/ - 62 + 6d)^ + (-6c- 661 - 18d + 362)- f 



^3 
Ttxi = 

T — 







4 2 

- 2061^ + (1261 + 262 - 2/ - 24d)^ 
t 



+ (66 - 4a - 461 + 362 - 24d + 2/)— + (126i - 66 + 48d 

X* 

^3 

+4a- 562)-^] 



Txx^ = 

T.,., =^[l06i^ + (/- 62 - 146i)^ 
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62 



5^8,1 , ^ 2^ 
3/+ 361)^ + (-4a + 62 + -/ 



-62 +46l)— ] 



/ , FIG. 4: (color online)The contours of energy density, in unit 

of as a function of t and xi. The magnitude is repre- 

sented by color, with darker color corresponding to greater 
magnitude. 

(38) 



We plot the energy density profiles in FiglH 

VIII. SUMMARY AND DISCUSSION 

The main objective of this paper was to calculate a 
"hologram" of the falling open string, which has ends 
attached to heavy quarks moving with constant veloc- 
ities ±u. After an appropriate tool - Green function 
for time-dependent linearized Einstein equation - was 
constructed, and stress tensor of the string calculated, 
a convolution of the two gave us the stress tensor of an 
"explosion" seen by an observer at the AdS boundary. 
Apart of analytical results in different limits, we have 
given pictures of the time evolution of the energy den- 
sity and the Poynting vector in FigslU [H In short, our 
main finding is that it looks like an explosion, with mat- 
ter "fireball" expanding from the collision point, but a 
non-hydrodynamical explosion, in which fluid cannot be 
assigned temperature or entropy. 

What can be a physical significance and applications 
of these results? 

Literally, they describe energy/momentum flow fol- 
lowing a collision of say two heavy-quark mesons in a 
strongly coupled Af—A gauge theory). It would then 
be instructive to compare these results with those in 
a weakly coupled regime of the same theory, in which 
the appropriate calculation would be perturbative ra- 
diation of massless gluons and scalars. Those are well 



known to produce dipole radiation at small velocities and 
bremmstrahlung cones (or "jets") in forward and back- 
ward directions. In QCD perturbative radiation is af- 
fected by confinement effects as well as the presence of 
light fermions: thus formation of QCD strings and their 
breaking by light quark pair production. All of it is well 
modeled by QCD "event generators" , one of which ~ the 
Lund model - we mentioned in the Introduction. 

So, why one would be interested in a strongly-coupled 
version of the "event generator"? One reason can be 
methodical: to better understand the difference between 
strongly coupled conformal regime and confining theo- 
ries, as far as jet physics is concerned. It has been studied 
in literature that hypothetical "hidden valleys theories" 
which may be found at LHC 19] are strongly coupled: so 
there is some nonzero (but tiny) chance that those can 
be used in real experiments one day. 

However, as explained in the Introduction, those were 
not our motivations. We have done this calculation as 
a methodical step toward understanding heavy ion col- 
lisions, in the AdS/CFT setting. One obviously needs 
to study one falling string before considering many. And 
simply adding the effect of many strings, as we did above, 
is not yet sufficient to understand heavy ion collisions. 

As a discussion item, we would like at the end of the 
paper to indicate where we will go from here. What 
we would like to understand in general is how and under 
which conditions the equilibration and entropy production 
happen, so that non-hydro explosion described above be- 
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comes hydro-like. In order to derive that, one has to 
abandon the "probe approximation" used above, and in- 
clude the gravitational impact of falling matter (strings 
in our setting) in the metric. Only then one may see a 
transition from extremal black hole (AdS metric we use) 
to non-extremal black hole with matter mass added and 
a nonzero horizon formed. The horizon, when present, 
does provide both Hawking temperature and Bekenstein 
entropy. We expect to use in the next paper of this series 
a "two- membrane paradigm", in which collision debries 
are represented by one (falling) membrane, and the (ris- 
ing and then falling because of stretching) horizon mem- 
brane by the other. (These two membranes can be asso- 
ciated with the so called top-down and down-up equili- 
bration scenarios, proposed in literature in various model 
settings.) When and where most of the entropy is pro- 
duced is the major issue to be addressed. How that is 
reflected in the "hologram" observed in the gauge theory 
could then be calculated in the linearized approximation, 
as above. 



sin{^x)xdx — TTYQ{a\/P — ^2) x 

Vx'^ + a'^ 



2 r ^^^^^^^^lsin{^x)xdx - -2K'Q{a^/W^) x 



^x'^ — a? 



- h) 



2 r '''^^^'^f^f\ in{^x)xdx - -iT4{a,/W~^) 
Jo + a 



.e{h-0 + '^6{b^i) 



(A3) 



APPENDIX A: INVERSE FOURIER 
TRANSFORM OF GREEN'S FUNCTION 

Let us recall the expression of Green's function in mo- 
mentum space: 

G(z, z') = { -frh{-iXz<)K2{-i\z >) w < 0, > k® 



jrh{\z^)K2{\z >) Ic^l < k® 



with a,b,S, > 

These formulas are obtained by differentiating with re- 



spect ^ the cosine-transform of 



and ^^^^iyf+f) 

Let's focus on the case t > at the moment. With the 
help of IA3[ [X2] becomes: 



(Al) 



We will use Q(2i(3) to refer to the three cases as in- 
dicated above. In order tojdo the inverse Fourier trans- 
form: / '-'(■^' z')e^^^~^^^ dw(i^k, we make the change 
of variable for each case: 



G{z, z')e"^'-"'''dujd-'k 



(2^) 



^ ^ G(z, z')e^^T^^^-f^edk^^ 



kr Vfc2 + A2 



1 



cj = -Vk^^——- I G{z,z')e"^'-"''=dLjd^k 
(27r)^ J 

1 



(27r)3 



Giz,z')e 



(27r)4 . 

,^„-tVmi:Tt '2sin{kr) XdX 

kr Vfc2 + A2 



= ±\ k^ - A2- 



1 



{2n) 



G{z,z')e"^'-'^^dujd^k 



i2n) 
XdX 

Vk^ - A2 



^ f G{z,z')2cos^Jk^-XH ^'''!^^''h ^dk 



kr 



(A2) 



Q: / -|/2(zAz<)i^2(«Az>)^ 



TTY^{X^/t^ - r^) 



=0(t-r) - 2i\:^(Av/r2 -i2) 



A 



Vr2 - t2 
8{t - r) 



0{r - t) 



e{t-r) + nr 

A rr>A f 



\Jt^ — r 



|_2|/2(-^Az<)i^,(-^Az>)^ 
=e{t~r)~2K'^{xVr^) 



TrY^{X^/t^ -r^) 
^ 0{r - t) 



.e{t -r)~in 



Vr2 - 
S{t - r) 



2 / -4/2(Az<)i^2(Az>)-^^^ 



(27r)3 



-A 



--0{t - r) + TvY^iXVr^ 



2K'„{XVt^ - r2) 
A 

-t)) 



Vr2 - i2 



(A4) 



We use the following formulas to evaluate the k- 
integrals: 



Suppose we replace z> by z> — ie, convergence of the 
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integral enables us to rotate the contour of (3) and (T) 



\^t^ — r'- 



-At-r) 



-—l2{iXz<:)K2{i\z>)--—— x 
z (Ztt)'^ 



If we are only interested in the coefficient of the 
termp2j. we may make the following substitution: 
l2i±iXz<) -lX^MXz<) ^ iA2,z> ^ z' 

Further evaluation of the integral involves the two for- 
XdX mulas: 



{-TTY^{Xy/t^ - r^) ~ mJ'^iX^Jt^ - r2)) 



e{t - r) 



-l;l2{lXz<)K2{lXz>)^^ 

X 



4/2(Az<)i^2(Az>)-^ X 
z' (ztt)-' 

A 



f^OO 

/ x''K^{ax)dx = 2i'-'^a-''-^r 
Jo 



(a2 + 62)^'+''+i 



(A8) 



(-^yo'(^Vr2 _ ^2) + i7rJ^(Av/r2 - _ t) 



Vr2 - t2 



It is easy to find / X^ dXK2(iiXzy) and 
/ A4dAif2(Az>)Ji(AVr2 - t^)9{r - i) contain no 



Similarly, suppose 2:> ^ z> + ie, we can rotate the 
contour of (3) and (2) 



(AS^ingularity in z>, thus the contribution from ± le 
cancel each other. We end up with: 



-^l2iXz<)K2iXzy)-^^ 



z' 
-A 



-2K'oiX\/t^ - r^) X 



-l^K2i^Xz^)-^{-2n^)MX^/?^) 



-At- 



= I --^/2(-«Az<)i^2(-jAz>)-^^ X 



A 



{~7TY^{Xy/t2 - r^) + iiTj'^iX^t^ - r2)) 



-A 



- r2 



(2^) 

0(t - r) 



x/i2~7: 



2> — '2:> —26 



2 , , , , AdA 

— l2{~lXz<)K2{-lXzy)j^ 



-2K[^{Xy/r^ - t^) x 



A 



AdA 
(2^ 
12iz' 
-J2^ 
x9{t-r) 



X 



2 -A^ 
z' 8 

e'(i- r) 



i^2(-iAz>) 



2:> — »-2:> +ie 



(^2 - r2 - z'2 + ie)4 (i2 - r2 - z'2 - ie)4 

(A9) 



(-^yo'(^v/r2 _ ^2) _ in^iX^r^ - t^)) 



^l2iXz<)K2iXzy)^^ X 

, For t < 0, similar procedure leads to a vanishing re- 
■ fl^T- — 3ult. Therefore the term of the Green's function in 
2 — f2 Jcoordinate space, which is exactl y th e propagator we are 

(A6)ooking for, can be expressed as [2a |: 



Summing up piece 0,(2) and (3) we find various terms 
cancel against each other. We are left with: 



Pr = 



z> ^ z> — le : 

-^l2{lXz<)K2{lXz>)-^ 

-A ^ d{t-r) 



12iz' 

1 



(t2 _ j-2 _ ^z: 



-2niJl^{X^/t^ - r2) 



(^2 _ r2 - z': 



6'(t - r) 



(AlO) 



AdA 



(27r)3 

z> ^ z> + ie 



^iJ'oiX^r^^) 



X 



Vr-2 - i2 
AdA 



+ / --/2(Az<)i^2(Az>) 



APPENDIX B: EXTRACT ENERGY DENSITY 
IN COMOVING FRAME 



-—l2{-iXz<:)K2{-iXzy, 
z' [It:)-^ 



The aim is to kill the (0, i) components(momentum 
density) by local Lorcntz transformation: T' — STS^, 
2Tri4{XVt^ - r2) wliere 



-A ^ S(t^r) 
X ft{f - r) - '- 



XdX 
(2^ 



+ / --/2(Az<)if2(Az>) 



^ = 7 ^ = ' Sati = (Bl) 



-■KiJ'oiX^Jr^ - i2) 



A 



Vr2 - t2 



Local Lorentz transformation matrix is a prod- 
uct of matrices, which are either rotations, e.g.Q or 
(A'/l)oosts,e.g.(2) This is a consequence of the fact: SgS^ = 
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where g = diag{—l, 1, 1, 1) is the Minkowski metric. 



cos(6l) sin(6') 
-sin(6l) cos(6i) 



1/ 



/chy shy 
shy chy 

\ 



(B2) 



A nice property of g is: g = = g ^, therefore gT' 
and Tg are related by similarity transformation: 



gT' ^gSTgg^S^ ^gS{Tg){gSy 



(B3) 



Since gT' does not alter the zero entries of T', now the 
problem becomes to kill the (0, i) components of Tg by 
similarity transformation. The original matrix Tg can be 
viewed as an operator L acting on a set of basis, while 



the similarity transformation is just a change of basis: 



(B4) 



{gT')mn = {e'^,Le'^) 

{gS^nm^in 

We want to find a basis c'q such that Lc'q 
A — {gT')oo. Denote Cq = XmSrm it is easy to show 
Xm{Tg)mn = Aa;„ This is exactly an eigenvalue problem 
for matrix (Tg)'^ . The restriction in transformation ma- 



Aeg with 



trix is translated to: 



> 0. It turns 



out this condition is just enough to determine a unique 
eigenvalue out of four possible eigenvalues. The energy 
density is given by e = — A. 
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the same boundary condition is used in [ij as a limiting 
case of thermal AdS background 

We have assumed that the pole lies in the integration 
range of t in using the residue theorem. In fact if it is not 
the case, we will still get a vanishing result due to the 
simple cancellation between ± term in P 
By analyzing the small z expansion of hmn{z) = 
/ G{z, z')smn{z')dz'd'^x, oue can show the coefficient of 
z^ term equals that of z% term, provided that Smn{z') 
contains no 2'" and z'^ terms. The latter condition is 
satisfied by the sources considered in this paper 
the t > condition is included in the theta function 



